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James Franklin Campbell 
ABSTRACT 

A preliminary investigation has been undertaken to theoretically 
determine the geometric and aerodynamic parameters which have the 
greatest influence on the stability and performance of the nonporous 
towed decelerator. To aid in this endeavor a mathematical model was 
generated to describe the flexible tow line - decelerator dynamical 
system as a rigid two-body problem. The resulting second-order gov- 
erning differential equations of motion were used to obtain the charac- 
teristic equation (quartic) describing the coupled motions of the two 
bodies. Evaluation of the coefficients of the quartic yielded expres- 
sions which illustrate when the i Tiamical system is unquestionably 
unstable . 

The two single-degree of freedom cases investigated were the 
motion of the decelerator about its nose and the motion of the rod 
and decelerator about the pivot point. For both of these cases 
expressions were derived for the natural angular frequency, damping 
factor, steady-state solution, and time to damp to one-half amplitude. 
Specific comments were made concerning a conical shaped decelerator to 
illustrate some of the geometric and aerodynamic parameters affecting 
the spring constant and the natural angular frequency. The boundary 
between stability and instability was obtained for these two single- 
degree of freedom cases. 
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TV. INTRODUCTION 


Advent of super- and hypersonic aircraft, missiles, and space- 
craft has dictated the need for adequate deceleration systems capable 
of performing emergency rescue and/or payload recovery. For example, 
rescue of a pilot ejecting from his aircraft or spacecraft at extremely 
high speeds is dependent on some decelerator system to stabilize him 
and to slow his fall to the point where a conventional terminal descent 
system, such as a parachute, can be utilized. Also, the success of 
a mission to place an instrument package on the surfaces of Venus and 
Mars largely depends on proper deceleration of the payload through 
their respective atmospheres. One such system being considered to 
perform these types of missions is the towed deployable decelerator 
which utilizes aerodynamic drag to decelerate and stabilize the payload 
or capsule. Two of the approaches made in designing towed decelerators 
required to function at supersonic and hypersonic speeds are: (l) 

adaption of existing subsonic parachute configurations and (2) genera- 
tion of nonporous inflatable bodies having specific geometric shapes. 

Of these two approaches, the nonporous body inflated either by 
compressed gas or by ram-air (inlets) presently appears to provide 
the more promising approach to high-speed deceleration problems, as 
well as having perhaps the greater growth potential for future 
applications . 

Some insight into the stability and performance of several 
nonporous decelerate'" configurations is provided In the experimental 
wind-tunnel investigations of references 1 and 2. The scale of the 
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test models, however, was severely limited by the size of the wind 
tunnel facilities, thus limiting the usefulness of the data. The 
degree of decelerator stability reported in these investigations 
was obtained by observation, a method typical of most wind tunnel 
tests of towed decelerator systems . Because of the inherent difficulty 
in obtaining adequate measurements of decelerator stability, emphasis 
has been placed by some investigators on the use of analytical 
techniques to determine decelerator stability (references 5 and 4) . 
Though the experimental investigations and analytical studies have 
illustrated some trends in decelerator dynamics, there still exists 
a general lack of understanding of the principal factors affecting 
the dynamic stability of a towed decelerator. 

In an effort to shed new light on this problem, a preliminary 
investigation has been made to determine the geometric and aerodynamic 
parameters which have the greatest influence on the stability and 
performance of the nonporous towed decelerator. One of the inherent 
difficulties arising in this type of study is the variable nature of 
the geometries of the tow cable and decelerator, both of which are 
flexible. Also, the towed decelerator is required to operate in the 
nonuniform flow field aft of the forebody which is being decelerated. 

As a first step in the analysis of this complex stability problem, 
several simplifying assumptions are made; the tow cable is assumed 
to be a rigid rod with a pin connection at each end and the decelerator 
is assumed to be a rigid body of revolution, the resulting dynamical 
system being described as a coupled two-body problem. Since the cone 
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is the basic shape about which other decelerator configurations are 
designed, specific comments are made concerning this shape during the 
course of this study. For the present analysis the stability and 
performance of the towed decelerator are examined for uniform (free 
stream) flow conditions > such conditions would exist when the wake 
of the forebody has essentially disappeared, corresponding to a very 
large tow cable length. Appropriate comments are made to relate the 
effects of the wake on the stability and performance of the 
decelerator . 



V. LIST OF SYMBOLS 


A 


*1 


B x' \ 


a 


*L 

5a 


C 


m 

dC 

m 

m ~ da 
a 


N 

c 

ckT 

a 

CG 

c 


aerodynamic axial force on decelerator 
reaction forces at attachment point between payload and 
decelerator system 

length to rod center of gravity, measured along rod axis 
from origin of rod 

length to decelerator center of gravity, measured along 
decelerator axis from nose of decelerator 
reaction forces at attachment point between rod and 
decelerator 
length of rod 

. , „ . , aerodynamic axial force 

axial- force coefficient, * = 

%> s 

drag- force coefficient, ^rodync-clc dreg force 
lift-force coefficient, ^odgamlc lift force 

lift- force coefficient slope at zero angle of attack 

, , , , , . , aerodynamic pitching moment 

pitching-moment coefficient, - ■_ — — 

q„s c 

pitching-moment coefficient slope at zero angle 
of attack 

normal- force coefficient, force 

normal - force coefficient slope at zero angle of attack 
center of gravity 

damping term in equivalent spring-mass-dashpot 
dynamical system 
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5 


c 


cr 


% 

D 2 

g 

h 

J 


k 


dL 

3a 

3l 

3d 



critical damping representing the limiting case 
between oscillatory and nonoscillatory motion 
reference length of decelerator, taken to be base 
diameter for conical shapes 
aerodynamic drag force on rod 
aerodynamic drag force on decelerator 
acceleration due to gravity 
height of conical decelerator 

sum cf mass moments of inertia of rod and decelerator 
about the respective centers of gravity of the rod 
and decelerator 

mass moment of inertia of rod about the center of 
gravity of the rod 

mass moment of inertia of decelerator about the center 
of gravity of the decelerator 
spring -onstant in equivalent spring-mas s-dashpot 
dynamical system 
aerodynamic lift force on rod 
aerodynamic life force on decelerator 

slope of lift-force curve at zero angle of attack 

rate of change of lift force with rate of change of 
angle of attack 

rate of change of lift force with pitch angle 
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oL 

dd 2 

*2 

3m 

3a 

3m 

3d 


3m 

3m 3m 


302 

m 


“l 

“2 

N 

3n 

5a 

q 

% 

Re 

Ro 

r 

and 

S 


rate of change of lift force with pitch velocity 

aerodynamic pitching moment on rod 
aerodynamic pitching moment on decelerator 

slope of pitching-moment curve at zero angle of attack 

rate of change of pitching moment wit v rate of change 
of angle of attack 

rate of change of pitching moment with pitch angle 
damping of pitch 

mass term in equivalent spring -mass -dashpot iynamical 
system 
mass of rod 
mass of decelerator 

aerodynamic normal force on decelerator 

slope of normal -force curve at zero angle of attack 
pitch velocity of decelerator, 
free- stream dynamic pressure 

Reynolds number based on decelerator base diameter 

Routh' s discriminant 

radius 

roots of governing differential equations of motion 
describing the transient (or complementary) solutions 
reference area of decelerator, taken to be the base 




area for conical shapes 
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t 

At 

V 

V. 


x, y 


V y i 

x 2 , y 2 


3 

7 

8 



time 

time to damp to one-half amplitude 
cone volume 

flight speed of payload- decel era tor system 
vertical velocity of decelerator center of gravity 
vertical velocity of nose of decelerator, resulting 
from rotation about the attachment point between 
payload and decelerator system 
vertical velocity of center of gravity of decelerator 
with respect to nose of decelerator, resulting 
from rotation about the attachment point between 
rod and decelerator 

cartesian coordinate system (X axis tangential to 
flight path) 

cartesian coordinates of rod center of gravity 
cartesian coordinates of decelerator center of gravity 
angle of attack of decelerator center of gravity 
angle between payload reference axis and horizontal 
angle between flight path and horizontal 
cone semiapex angle 

angle between rod axis and tangent tu flight path 
angle between decelerator axis and tangent to flight 
path 

induced angle of attack resulting from vertical 
velocity of the decelerator 1 s center of gravity 
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p 




CD 

n 


exp 


particular (steady state) solution of the differential 
governing equation of motion for decelerator rotating 
edjout its nose 

damping factor (ratio of actual damping to critical 
damping) 
density 

damped angular frequency 
natural angular frequency 
exponential 


In 


natural logarithm 


z 

8 

. (dot) 

. . (double dot) 


summation 

partial derivative 

first derivative with respect to time 
second derivative with respect to time 


Subscripts 

1 induced angle 

Z large 

0 initial conditions (time equal zero) 

p particular solution 


s 


small 


VI. THEORETICAL DESCRIPTION OF TOWED DECELERATOR 


Mathematical Model 

This analysis was performed in order to determine the geometric 
and aerodynamic parameters having the greatest influence on the 
stability and performance of the towed decelerator system. In the 
most general case the dynamics of the towed decelerator are necessarily 
coupled with the dynamics of the payload, and vice versa, making a 
straight forward analysis of the decelerator * s dynamic characteristics 
rather difficult. In this investigation, therefore, the angular 
motion of the payload is assumed zero, so that the problem is sim- 
plified and the following analysis can be concentrated on the towed- 
decelerator dynamics. Since orientation of the flight path does not 
affect the stability of the decelerator, the decelerator system will 
be considered in a horizontal 
mode of operation (7 = 0 in 
Figure 1 ) . It should be noted 
that the most common method of 
experimentally testing towed 
decelerators utilizes a payload 
mounted rigidly in a horizontal 
position. 

In reality, the tow line and decelerator are both flexible and 
represent a multidegree of freedom problem. A general approach to 
describing this system would be to consider N number of particles 



Figure 1.- Schematic of payload- 
decelerator system. 


9 


10 


having individual masses. The more restrictive approach of assuming 
the tow line and decelerator as rigid bodies of revolution is utilized 
in the present analysis. Only plane motion is considered, yawing 
and rolling moments being neglected. 



Figure 2.- Two-body system. 

The towed decelerator system is illustrated in Figure 2, where 
pin A represents a fixed pivot point and has no motion, and where pin 
B is free to rotate with the system. Attachment of the rigid rod 
and decelerator is at the nose of the decelerator (pin B) . This 
resultant two-body problem is strikingly similar to the double- 
compound pendulum described by Ifyklestad in reference 5 . The basic 
differences between the two dynamical systems are: (l) weights in 

the present analysis represent near constant forcing functions for 
the system, whereas they act as spring constants in the double- 
compound pendulum problem, and (2) aerodynamic forces and moment 
are included in the present analysis. The method of obtaining the 
differential equations of motion for the present study is identical 
to that performed by Myklestad. The cartesian coordinate system 
employed for the two bodies is illustrated in Figure 3, and the 
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resulting free-body diagrams in Figures 4 and 5 . The X-axis is 
oriented to coincide with the direction of the free- stream velocity 
vector. The aerodynamic forces and moments of the rod and decelerator 
are shown in the figures to act at the respective centers of gravity. 



Figure 5*- Free body diagram of decelerator. 


In the analysis to follow, the aerodynamic forces and moment on the 
rod are assumed to be negligible, and likewise the rod's aerodynamic 
damping . 


Equations of Motion 

First, the internal reactions and are found from summing 

forces in the X and Y directions to be. 


^ + Dig Xg - Dg = 0 

Bx = D 2 " m 2 X 2 


( 1 ) 


and 
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By + m 2 y + L 2 - m 2 g = 0 

( 2 ) 

** = - l 2 - ”2 y 2 + ”26 


Two equations of equilibrium can be obtained from the free-body 
diagrams by taking moments about pins A and B. Summing moments about 
A yields 


^ M A = 0 = n^ga.^ cos ^ + m^h^a^^ sin 0^ - rryy.^ cos 9 1 


J 1 9 1 


- W sin 9 1 + Vi 


cos 9 


Substitute for B x and EL. 


from equations (l) and (2) to get 


X M A = C cos 0 2 * sin 9 1 “ m l^i a l cos 9 1 " J 1 9 1 

- b^Dg sin 0 1 + b 1 sin ^ m^ - b^ cos 9 1 m^y^ 

+ b^ cos 9 1 L 2 + n^gj (5) 

Summing moments about point B yields 

^ Mg = 0 = ra 2 g& 2 cos 0 2 + sin 9 2 “ *2*2*2 C0S ®2 " L 2 a 2 COS 9 2 


- D 2 a 2 sin 9 2 - J 2 0 2 + Mg 


( 4 ) 
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It is desirable to change the cartesian (CG) coordinates to the 
generalized coordinates of 9^ and 9^. To accomplish this the 
following relations are used: 

x x = cos 9^ x 2 = b^ cos 0^ + a 2 cos 9 2 

y x - sin 9 1 y 2 = b 1 sin d ± + a £ sin 0g 

Differentiating with respect to time yields 

*1 ■ - sin 9 1 x 2 = - b^ sin 9 1 - a 2 § 2 sin 9 2 

y x = a^ cos 9 1 y 2 > b 1 0 ]L cos 0 1 + a^ cos 9 £ 

= - a^S^ sin 9^ x 2 = - b.^ sin 9^ - b^B^ 2 cos 9^ 

- a^ 2 cos 9 ! - * 2*^2 Sln 9 2 “ a 2 ^ 2 2 C0S d 2 

'^1 = ^^l C0B 9 1 yg = COS 9 1 " ^l^ 2 S ^ n 9 1 

- a^ 2 sin e 1 + a 2 9 2 cos 9 2 - a 2 0 2 2 sin 9 g 

Substitute the appropriate expressions into equation "3) to get 
nijga^ cos 0^ + n^a^ sin 9^ (- a^9^ sin 0^ - a^0^ cos 9 1 j 

- sy^ cos 0 1 (y'y cos S x - y^ 2 sin 9 1 ) - 


( 5 ) 
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- b^D 2 sin 0^ + sin 9^ sin 9^ - b.^ 2 cos 

- a 2 9 2 sin 9 2 - a^ 2 coa 0 2 ) - cos (b^ cos 0 1 

- b 1 0 1 2 s in 9 ! + fi^ 9 2 cos 0 2 " a 2®2 2 Sin 9 2) 

+ b ± cos e i (- L 2 + n^g) = 0 
Assume small displacements so that 


sin 9 *» 9 and cos 9 *» 1 


• 2 *2 

Higher order and nonlinear terms such as 9^ , and 0 2 are 

neglected; that is, products of the angles and products of their 
derivatives. Thus, 


- w + vi ( a A) + J A + va + i i m 2 (Vi 


- \ (- L 2 + “ 2 ®) = 0 


Collecting terms, 

2 \ ** 

+ m 2 b 1 + JJ 0 1 + \^2 d l + a 2 b l m 2 9 2 + b l L 2 = “l 8 ®! + m 2 gb l 

( 6 ) 

Substitute expressions (5) Into equation (4) to get 

. . O 

m 2 ga 2 C0S ®2 + “2*2 Bln 9 2 b l 0 l sin 0 l " b l 0 i cos 9 i " a 2 0 2 sin 0 2 
- a^dg 2 cos 9g) - mga 2 cos 9 2 ^b.^ cos 9 1 - b.^ 2 sin 6^^ 
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+ cos 0 2 - a 2 9 ? 2 sin 0 2 ) - L 2 a 2 cos 0 2 " D 2 a 2 sin 0 2 

- J 2 0 2 + Mg « 0 

Again, assume small displacements (sin 0** e and cos 0 ** l) and 
neglect nonlinear terms. There results 


" m 2 ga 2 + ®2 a 2 ( b A + a 2^) + L 2 a 2 + D 2 a 2 9 2 + J 2^2 ‘ M 2 = 0 


Rearranging, 


m 2 a 2 b l 9 l + ( In 2 a 2 2 + ^ 2 ) ®2 + D 2 a 2 9 2 + L 2 & 2 " **2 = m 2 ga 2 (?) 


It is now desirable to consider the aerodynamic forces and moment 
existing on the decelerator. The reader should note that the subscript 
2 associated with the decelerator* s aerodynamics will be deleted during 
the remainder of the discussion. 

The angle of attack of the decelerator* s center of gravity is 
defined as follows: 


where 9 ^ represents the angle between the decelerator axis and the 
X-axis and 9 represents the induced angle of attack resulting 
form vertical velocity of the CG of the decelerator. 

Note; The vertical velocity of the decelerator CG consists of 
two components, v^ and v 2 , resulting from rotation 
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about pin A and B, respectively (see Figure 6). The 
component is the velocity of the decelerator CG with 

respect to pin B. 



Figure 6.- Induced velocities. 


Since small angles of attack are being considered, it follows that 

D , v __ v V 1 + V 2 

8 tan y ~ y _ 

x 00 00 00 


It caii be seen from Figure 3 that 


V 1 * Vl and v 2 = a 2 0 2 


Thus, 


9 _ 

2 i = V„ 


The angle of attack of the decelerator is then given by 


„ . Vi . “a e 2 

a - e 2 + — + ~T~ 


00 00 


( 8 ) 


Also, 


A , b l 8 l a 2 0 2 
“ = 9 2 + — + — 
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For this analysis the aerodynamic drag force (D) is considered to be 

2 

independent of angle of attack since D ~ a . The lift force (L) 
can be expressed as 


L = 


OL 

\ CL 

da 


8l . 

5d a 


In keeping with the assumptions of linearization, higher-order 
terms are not considered. Substituting a and a into the 
expression for lift force. 


8l 9 b l 8 l a . a 2 8l • 8 l a .1 8l •• a 2 8 l " 

L = 5E 9 2 + — SS 9 i + r^ 9 2 + 53[ 9 2 + v“^ e i + v"SdL 9 2 

00 OO 00 ULL 00 


( 9 ) 


The aerodynamic pitching moment can be expressed as 


M 8m 8m . 

M - -r — <X + T. T - yr— 

5* 8<i 5C 


8m a 8 m 

°p + —r~ 

2 d 89 2 


Substituting for a and a. 


dM 


b l 8 m 


M = 5^ 9 2 + iL^ 9 i + r 


00 


2 OM A . dM ’ , b l dM •• a 2 dK " dM * 

3 S 2 + as 9 2 - \r 5 - 9 1 + r 34 0 2 + d F 2 


( 10 ) 


where 


8m 8m 



Substitute the expression for lift force into equation (6) . 
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(V^ + "2V' + J l) + b l M l + a 2 b l“a 0 2 + b ! i 9 2 + Y- I 0 


a 2 b l bh • , b l bL 'A . a 2 b l dL '* 

+ — ss e 2 + b l s 5 2 + V7 55 9 1 + — M 2 = "I s 1 


+ nu<rb 




Collecting terms, 


("l 3 ! 2 + "2 b l 2 + J 1 + £ |) *L + Jr I S 1 + b l M 1 + (^Vs 

. a 2 b l bL\ *• . ( a 2 b l bL . . dL \ • dL p 

+ ^ ^ + l— ^ + b l ^) Q 2 +b l^V V*L 


m 2 gt 1 


for lift force and pitching moment into 


_ b 

2 a 2 b l 9 l + fa>*2 + J 2) 9 2 + Da 2®2 + a 2 I 9 2 + HP I! 9 1 

OO 

o o 


♦ |fc 5 + , to 5 ♦ *4± ^ ? + !S- -4 - i* 9, 

V. 3S 2 2 a£ 2 V TO 1 V* a<x 2 ^ 2 

b l dM • a 2 dM • SM; b l SM v a 2 dM v 

9 i-^3H e 2-ai 9 2-7:3l 9 i-^aZ ■ 


1 bM • 2 dM * dM 

’ e i " V M 3£ 9 2 " di 


dM * 

s; e 2 - ”2« a 2 


Collecting terms. 
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+ &gbl 
V. dA 


_i Sm 


9 1 + 


a 2 b l SL b l 5 m\a ( Z 

X" 5S - ^ 5S/i + (Va 


+ J, 


+ 


2 

a 2 dL 



S 2 cimVa ^T a 2 , SL a 2 dM SM dM\- 

2 l v^^ + a 2^“ 2 


+ ( a 2 D + «2 si - si) 9 2 ' “2 ea 2 (12) 

Thus it is seen that equations (ll) and (12) represent a system of 
second order, nonhomogeneous , linear differential equations with 
constant coefficients. 

Before proceeding to the discussion of the possible solutions 
to these equations, a few comments are necessary concerning the aero- 
dynamic forces and moment assumed to act on the decelerator. As 
seen in equations (9) and (10) not all of the factors influencing 
lift force and pitching moment are taken into account. It is 
believed, however, that the exclusion of these terms will not have 
an appreciable effect on the results of this analysis. It should be 
emphasized that the equations of motion apply for any decelerator 
shape as long as the shape is a body of revolution. One such body of 
revolution often considered in decelerator problems is the conical 
shape (this shape forms the basis for most decelerator designs). 
Experimental data obtained on a series of cone bodies at free- stream 
conditions (reference 6) are shown in Figure 7 and illustrate the 
linear nature of lift force and pitching moment with angle of attack; 
also, drag force is approximately constant at angles of attack near 
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2 

zero since D ~ a . Such conditions would exist on a towed conical 
decelerator when the wake of the forebody has essentially disappeared, 
corresponding to a very large tow cable length ^b^J . The question 
arises, however, of how the aerodynamic forces and moment vary with 
angle of attack if the cone is subjected to the nonuniform flow field 
associated with the wake. The answer to this question appears 
formidable at the present time because of the dependence of the aero- 
dynamics on the size of the cone (with respect to the forebody), 
the location of the cone in the wake, forebody geometry, free-stream 
Mach number and Reynolds number, boundary-layer buildup and separation 
effects on the tow cable, and on cone semiapex angle. The possibility 
exists that to adequately describe the aerodynamics of the cone in a 
nonuniform flow field, it may be necessary to retain higher order 
terms, thus defeating the simplicity of the approach of the present 
analysis. For these reasons, free-stream flow conditions are assumed 
to exist. The reader should note that the following analysis, though 
concerned principally with the conical shaped decelerator, could be 
applied equally well to any shape body for which '.he geometric and 
aerodynamic parameters are known. 
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shown in the stability investigation of reference 7 to be a justifiable 
assumption at high Mach numbers. Equations (ll) and (12) therefore 
reduce to 


(v 1 / + + j i) \ + v^- i + b i D9 i + (v*) 8 


f 21 a . . c>L 


9 2 + b l 5£ 9 2 = V a i + m 2 gb 1 


(13) 


and 


f . \ ^ ( a 2°l dL \ bu\ ‘ / 2 _ \ / a 2 5L 

( a 2 b l m 2) 9 1 + r — 35! - T 3S,J -1 + (®2 a 2 + J 2j 9 2 + Vv“ da 

a 2 dM ^ f t, ^ dL q 

" \T 3£ - 35 ) § 2 + + a 2 SS “ 3s] 9 2 = V a 2 (l4) 


Characteristic equation .- To arrive at a characteristic equation 
representing the system described by equations ( 13 ) and (l4) it is 
beneficial to write these equations as 


A 0! + B § + C 9 1 + D 9 g + E 0 2 + F 0 2 = G 


H 9 1 + T 0 1 + J 0 2 + K 0 2 + L 9 2 = M 


(15) 


where 


2 2 

A =• m^a^ + “g 13 ! + 


r> a ’ 2 ° 1 iL 

E = VT3S 


B _1 dL 

b -v;3h 

C - b x D 


D = a rt b 




F = b 


dL 
1 5a 


G = 


n^ga^ + n^gK 
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II = a g b m,. 


K = 


2 dL ^2 dM dM 
V* Sa " y x ckL~'Sq. 


T = 


a 2 b l dL _ \ dM 
V 3ci V 3a 


L = a,. 


(Is + D ) ■ t£ 


J - Va + J a 


M * "\g a 2 


Consider the homogeneous case. 


A 9 i + B 0 i + C 0 i + D V E 0 2 + F ®2 = 0 


h 0i + 1 §i + j e 2 + k e 2 + l e 2 = ° 


Assume solutions of the form 


ft O A a o At 

- e j e and 0 g = 0 g e 


so that 


0 2 = X 0! e 


At 


0 2 * A 0 2 e 


At 


01 = A 2 0i e At 


\ = A 2 9 2 e At 


Substitute and divide by e^ b to get 


A A 2 ©i + B A 0 X + C ?i + D A 2 0 2 + E A 0 2 + F 0 g = 0 


H A 2 0i + I A ?i + J A 2 0 g + K A 0 g + L 0 = 0 


(16) 


Collecting terms. 
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(aa 2 + ba + c) 3 + (da 2 + ea + f) o 2 = o 

(ha 2 + Ia) ? + (jA 2 + KA + l) 0 2 =0 

These are algebraic, homogeneous, linear equations for the constants 
9^ and 9 . They always have at least the trivial soltuion 

fy 

9 1 = 9 2 = which corresponds to = 0. These equations have 

a nontrivial solution precisely when the determinant of the 
coefficients is zero. 

(aA 2 + BA + c) (dA 2 + EA + f) 

(HA 2 + IA) (jA 2 + KA + l) 

Expanding this determinant and collecting terms yields the 
characteristic equation, 

(AJ - HD) A^ + (AK + BJ - HE - ID) A 5 + (AL + BK 4CJ - HF - IE) A 2 
+ (BL + CK - IF) A + CL = 0 (17) 

Numerical constants could be introduced and exact roots found from 
this equation, but for the purpose of defining stability it is 
necessary to look only at the signs which these constants carry. 

Stability criteria .- The necessary and sufficient condition 
for the stability of the dynamic system represented by t.ia quartic, 
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A' X + B' 


X 5 + C’ 


2 

X + D' 


X + E' =0 


is that 


and 


A', B' , C', D', and E' > 0 


R = B'C'D' - B ,2 e' - A'D’ 2 > 0 


This latter grouping of constants is known as Routh's discriminant. 
It is first necessary to evaluate the coefficients of the quartic, 
equation (17) > and to see under what conditions they might change 
signs . 

A' = JU - HE . (el^ 2 + ft,^ 2 + Jj] (ny./ + Jg) - 

= m 1 m 2 a ^ 2 + m 2 a 2 2j i + J 2 (^ m l a l 2 + m 2 b l 2 + J i) 


Thus, this coefficient is positive for all combinations of masses, 
mass moments of inertia, tow line lengths, and center of gravity 
locations. 


B' s AK + BJ - HE - ID « ^ny^ 2 + m^ 2 + jJ 


dL 

5a 


_2 oiA 5 m 
V Sa ’ 3q 


D 

b l 5L /_ . 2 . T \ / a 2 b l ( & 2 h l 5l \ 5m\ b _ 

V„ 55 ( m 2 a 2 J 2j ■ a 2 b i m 2 v, 55j ■ ^ V, 55 ' v w 55j a 2 b i m 2 


Expanding, 
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( 2 2 2 _ . 2\ 1 dL f 2 x T \1 

" ' K a l a 2 + J l a 2 + J 2 b l ) r 35 - (Vl a 2 + J l a 2 - 35 
\ '00 ' / 00 


- (vi 2 * “A 2 + J i) i 


^ and 
oa 

dM 

dq are 

bL > n 
da ~~ 

When 

and is given by 


the expression 


(VlV * J l a 2 2 + J 2 b l 2 ) sl * (VlS + j l a 2) 3S 

+ V » (“l 8 ! 2 + Vl 2 + J l) 3i 


(18) 


Also, 


+ h* 81 

+ v“~dd 


00 


a 2 dL a 2 dM dM 

V 3a’v 3a dq 

QCj on * 


\ r (*l . ' 

\ dM 

J |^2 + D 

j ' da 


+ + J 2 ) • a 2 b lS 35 


a 2 b l dL / a 2 b l dL ^1 dM^ 
V oo 35 v V » 35 " v oo ^ 


Expanding this expression and at the same time noting that ^ ^ + D 

(see subsequent discussion on page 42), it can be shown that 

C = m 1 a 1 2 a 2 + m 2 a 2 b 1 2 D + a^ || + b^^a^ 2 + Jg) - g| 3~ 

’ H (“l^ 2 + * 2*1 + J l) 


29 


wit 

As seen in Figure 7 for conical bodies, ^ is always greater than 
zero. Thus, when ^ > 0, this coefficient is definitely positive. 
V.'.en ^ < 0, C* is probably positive since the negative term 
involving ^ would be outweighted by the other positive terms. 


D 1 = BL + CK - IF = r 


1 oL 


V 3a 


(®s (is + D ) ■ si + ^ 


_2_ dL 
V 5a 


^2 3m _ 3m\ 3l [ a 2 b i 3 l \ 5m 
5a " 5qf 1 5a l V 5a ’ V 5E 


which expands to give 

r,. !. V 2 J. . K. D T> Sm v r, 

d . ^ ^ 55 - T- B ss - V 35 

x '00 00 

When the lift-curve slope is greater than zero, the coefficient D’ 
is definitely greater than zero. When D* = 0, the lift-curve slope 
is of negative value and is expressed as 


5a 7a~+TT 5a a 2 ( & 2 + b^ 5q 


Lastly, 


E’ = CL = 


= b i D a 2 (si + D ) " 


v n f a 3n 3m' 

= i D {*2 SS ■ 5S, 


Thus, this coefficient is always positive. 
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It has been mentioned that it is a necessary condition for the 
coefficients of the stability quartic to be positive valued to ensure 
a stable dynamical system. If any of the coefficients, say B' and/or 
I)', were to become negative, then the dynamical system is unques- 
tionably unstable. This has been shown to be the case when the lift- 
curve slope has a negative value less than that expressed by equations 
(l8) and (19) . Some support for this argument in the form of 
experimental data has been published in reference 1. These data 
indicate that towed cones with semiapex angles up to 4o° are definitely 
stabie, 45° cones being marginly stable, and cones with larger semi- 
apex angles being unstable. The major change in any of the aerodynamic 
parameters between a 40° cone and a 50° cone is the change in sign 
of the lift curve slope (see Figure J), being > 0 for the 40° 
cone and ^ being < 0 for the 50° cone. 

To show definitely that the quartic, expressed in equation (17), 
represents a stable system it is necessary to analyze Routh's 
discriminant, 

R - B'C'D f - B' 2 E' - A'D' 2 
o 

However, manipulation of the coefficients to provide a general 
expression for Routh's discriminant without making any simplylng 
assumptions or without substitution of numerical values is extremely 
tedious. Though evaluation of this discriminant is left to future 
investigations, it is believed that Routh's discriminant will be 
greater than zero if the lift- curve slope is greater than zero. 
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Descartes* Rule of Signs .- To determine the number of real roots 

existing for the quartic, equation (17), Descartes' Rule of Signs is 

used and is stated as follows: 

A polynomial equation, f(x) = 0, with real coefficients and 
arranged in descending powers of x, can have no more positive 
real roots than there are variations of sign between successive 
terms in f(x), and can have no more negative real roots than 
there are variations of sign between successive terms in f(-x). 

(reference 8) 

As has been shown, the quartic of equation (17) may be wiitten as 


l'A 1 + B'A 5 + C'A 2 + D* A + E* = 


where A', B', C', D', and E* are real coefficients and are equal 
to those constants shown in equation (17). If 0, it has been 
shown that all the coefficients in equation (20 ) are positive. 


A* ,B' ,C* ,D* ,E* > 0 

Applying Descartes' Rule of Signs it is seen that for f(A) all the 
coefficients are positive so that there is no variation of sign. Thus, 
there can be no positive real roots. For f(-A), equation (20) becomes 

A'A 4 - B'A 3 + C’A 2 - D'A + E* = 0 

showing that there are four (k) variations of sign. Therefore, there 
can be up to four negative real roots. Since there are no positive real 
roots, there will be either zero, two, or four negative real roots. 

If ^ < 0, it has been shown previously that A',C',E' > 0 and 
3' and D* can be less than zero. Letting B' = - B* ' and 
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D‘ = - D' ' in equation (20), 

A'A^ - B"A ? + C’A 2 - D"A + E’ =0 (21) 

Applying Descartes’ Rule of Signs it is seen that for f(A) there 
are four (4) variations of sign indicating that the quartic can have 
no more than four positive real roots. For f(-A), equation (2l) 
becomes 

A'A 4 + B' r A 5 + C'A 2 + D* 'A + E' =0 

showing that there are no variation of sign. Therefore, there can 
be no negative real, roots. Thus, if any real roots er.ist, they will 
be positive indicating instability. 

Review of Solutions to Second-Order Differential 
Equations and Stability Criteria 
It is now desirable to consider two particular cases of single 
degree of freedom. These two special cases are (l) motion of the 
decelerator about its nose (pin B) with 9^ = 0°, and (2) motion of 
the rod and decelerator about the pivot point (pin A) with 0^ = 0^. 
Since both of these special cases are represented by a second-order 
differential equation, it may be beneficial for the reader to review 
the types of transient motions possible with this type of equation. 

The equivalent spring-mass-dashpot dynamical system is described 
by the equation of motion, 

m 5 + c 0 + k 0 = 0 


( 22 ) 
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the roots of which are given by 

o „ c 2 k , „ , r l* , „ r 2 t 

r l,2 ‘ * 55 1 ' m uhere 9 - *1 e + *2 6 

The limiting case between oscillatory and nonoscillatory motions is 

when the radical equals zero. The value of damping for this condition 

is called critical damping, c . 

Cl* 


= 4 mk or c = 2 Jmk 

r*Y* * 


and since 


u) = .1 — , then c = 2m cd 
n \ m * cr n 


The dimensionless ratio c Ic is called the 

I cr 


bor, | . 


c c actual damping 

~ oT * critical damping 


cr 2 >/im “n crltlcal damping 
Using these definitions, the original equation (22) can be written 


in the form 


9 + 26 <d 0 + cd 0 = 0 
a n n 


where the root '; are now given by 


r l,2 " “n ± “n V - 1 


It is seen that the roots are: 




(a) real If § > 1; the motion is overdamped and the transient 

solution is defined by two exponential terms 

(b) equal if 5=1} the actual damping is equal to the critical 

damping 

(c) complex if 6 < 1} the motion is underdamped and the transient 

solution is oscillatory 

(d) imaginary if 5=0} oscillates at natural frequency with no 

damping. 

As a consequence, 5=1 separates the regions of motion (see Figure 8). 



This figure also illustrates one of the basic characteristics of a 
dynamical system: the motion can be classed as stable provided the 

coefficients of the governing equation of motion (23) are greater than 
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zero. Thus, 

26 a) =->0 and o > 2 = - > o 
n m n m 

This simply states that the damping and spring constants have to he 
positive valued to insure a stable system. For a more detailed review 
of the transient solutions of second-order differential equations 
and the resulting stability criteria, the reader is referred to 
references 5 and 8. 

Now that the types of transient motions possible for a dynamical 
system have been reviewed, the special case of the motion of the 
decelerator about its nose (pin B) will be considered. 



* niggag 


Assume a solution to match the right side of this equation. 


Therefore, 


0 2 = constant 

P 


S 2 - S 8 .0 

P P 


( 1 § + D ) “ si] 9 2 = m 2 ga 2 
v J P 


The particular solution for equation (2-), giving the equilibrium 


position of the cone, is thus 


m 2 ga 2 


v fe (§♦»)- £1 

The homogeneous condition of equation (24) is seen to be 
/ 2, t \ a xf^ ^ a 2 dM SlAi . fb L X J\ dM 

( m 2 a 2 + J 2) 0 2 + \J£" ’ SqJ 0 2 + |^2 ^35 + D J “ Sa_ 


This is in the same form as equations (22) and (23) in the review. 


where 


(tl + D ) - i! 


"a* i *\ ST — 

\ mgSg + J 2 


I = 2 — = 2 
v/km. 


a 2 SL ^2 SM SM 
V 5a ' V 5a Sq 


.a 2 + J J fa„ + D) 



37 


As pointed out in the review of solutions to second-order differential 
equations, the types of roots depend on the value of the damping factor, 
4* The transient (or complementary) solutions associated with the 
homogeneous equation, (26), for various values of damping factor are 


derived in the following discussion and are added to the particular 
solution, equation (25), to obtain the complete solution of equation 

(2k). 


> 1 «- The transient solution of equation (26) for 4 > 1 is 


r i t r 2 t 

S 2 = A e 1 + B e 


where 


r 


1 


-4 +co 
n n 



and 




The transient solution can be written as 

9g = A exp (-5 <u n + \j{ 2 - 1) t + B exp f-S - l) t 


e 2 = « 


-4o> t 


A exp 1 a> n t) + B exp J ^ - 1 o^t) (27) 


Combining equations (25) and (27), the complete solution of equation 
(2l+) for 4 > 1 is seen to be 


9 2 “ e 


n A exp (^4 2 - 1 o> n t) 


+ B exp V5 " 1 “nV + 0 2 


where the constants A and B are determined by initial conditions. 


The motion of 9 ^ as a function of the dimensionless parameter 
u> n t is shown in Figure 9 for several values of damping factor. As 
would be expected, the motion is an exponential function of time, the 
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more sluggish motion occuring for the larger damping factor. To 
obtain the solution for 5=1 it is necessary to again utilize 
equation (26). 

6 = 1 .- The roots of the characteristic equation are 


so that the solution is 


-CD t -CD t 

0£ = A e n + Bt e n (29) 

Combining this with equation (25), the complete solution of equation 
(2^) for |=1 is 

-CD t 

0 2 = (A h Bt) e n + 9 g (30) 


Note that 9 2 ->0 2 as ^ “* eo * T* 11 - 8 is the expression used to 

P 

illustrate the motion of 9^ with cu^t in Figure 9 for 5 = 1. 

6 < I .- The type of dynamical motion occuring more often than 
not in aerodynamic stajility problems is that which is oscillatory 
in nature, having a value of actual damping less than the critical 
damping (5 < 1) . In this case the roots of the characteristic 
equation are complex and are written as 

r l = “5 + i^ n Jl - S 2 and r 2 = -5 cD n - ico n \jl - % 2 

The solution of the homogeneous equation (26) is thus, 

0 2 = e ^ n (a cos \|l - i 2 - o> n t + B sin Jl - i 2 «Dt j (31) 

2 

where cd^ = cd^ ^Jl - { is the damped frequency of the motion. 





4o 


Combining this transient solution with the particular solution, 
equation (25), yields the complete solution of equation (24) for 
I < 1. 



The motion resulting from this expression is shown in Figure 10 for 
several values of damping factor. The largest value of 5 results in 
the motion being damped at the least value of o^t; the smaller the 
damping factor the larger the amplitude at any given time and the 
longer it takes the motion to damp. This trend is more easily noted 
by examining the time to damp to one-half amplitude. 

Time to damp to one-half amplitude .- Using the exponential part 
of e uation (32), the exponential envelope of the oscillatory motions 
seen in Figure 10 is obtained as 




where C is a constant. Letting 9^ correspond to a time t^, and 
correspond to a time t^, it is desirable to determine the time 


necessary to damp to a condition where 



Thus, 


V 1 

e 


4 1 


Define At = t 0 - t^ as the time to damp from an amplitude of 9^ to 
an amplitude of O^'. 



At 


and 


At = 


_ In l/2 _ - In 2 _ In 2 




-*«> 




n 


Recalling that 


a 2 5L ^2 <^M c)M 

V 5a V 5a " Sq 


*“n = 5 = 


<V2 2 J 2 ) 


the time to damp to one-half amplitude then becomes 

2 


At = 


2^a 2 2 + In 2 


a 2 5l _2 dM 5 m 
V 5a V 5a 5q 

00 00 


(33) 


(34) 


Spring constant .- It has been shovn that a criterion for 
stability is that the spring constant be greater than zero. From 
the governing equation of motion for the decelsrator rotating about 
its nose, equation (24), the spring constant is seen to be 



This expression can be rewritten as 


k = (^S c 



^4 




2 \ 


c 

=- C T + c_i 

- C 


? 1 L D 

m 


L. v. a , / 

a_ 


(35) 


This can be simplified by considering the lift force in terms of 


normal and axial forces. 


L = N cos a - A sin a 


dL 3N „ . dA . 

■t— = cos a - N sin a - t— sin a - A cos a 
da da da 


Recalling that drag was assumed to be constant with angle of attack 
(D ~ a 2 ), and for small angles drag force is equivalent to axial force, 
then = C. The above expression is evaluated at a = 0 ° to get 


dL dN 

3 a“S5.- A 


da da A 


where is equivalent to a* a = 0 °. Substituting this into 

equation (35 ) > the spring constant is found to be 


a 2 

■r- CL. - C 
c N m 

a a J 


To obtain some feel for the effects of the geometric and aero- 
dynamic parameters on the spring constant, equation (36) will be 
examined considering a conical shaped ch c^lorator. It is desirable 
to determine the difference in the spring constant for two different 
sizes of cone, the cone semiapex angle being a constant. It is noted 
that, since the semiapex angle is the same for the two cones, then 




the aerodynamic characteristics, and , are the same. Also, 

a a 

the value of a n /c would he the same. 


s l 


s = small 


l 3 large 


Writing the spring constant for both sizes of cone. 


k = a S c — C_ - C 
s ^s s s\ - N m 
Vc a a 


k. = a, s, c, — c„ - c 

l l l l\ - N m 
\c a a 


Thus, 


a S c 
T3 s £ 


q ! S I °] 


<1, S, 0 

k » - rrr ^ (57) 

s s s 

Therefore, knowing the aerodynamic spring constant of a cone having a 
given size and semiapex angle, the spring constant for a cone of 
corresponding semiapex angle hut different diameter and/or dynamic 
pressure can be determined. For example, consider the larger cone 
with a diameter twice that of the small cone (i.e. c^ = £ c s ). 
Assuming the same dynamic pressure (q ^ = q^). 


it is seen that 


u 


where 


and 


Therefore, 

k = 8 k 
Z s 

Thus, doubling the cone diameter for a constant dynamic pressure 
increases the aerodynamic spring constant by a factor cf eight (8) . 
If a particular size cone were considered in equation (37 ) > it is 
seen that doubling the dynamic pressure doubles the aerodynamic 
spring constant. This is important since the wake has a lower 
dynamic pressure than the freestream. If the cone aerodynamics are 
assumed unchanged from freestream to wake, then equation (37) 
predicts a lower spring constant in the wake than in freestream. 
Another factor to be considered concerning the spring constant is 
that it will remain the same magnitude regardless of center of 
gravity location. This can be seen from equation (3 8 ), where C 

m a 

is the aerodynamic pitching moment referenced to the decelerator 
center of gravity. To transfer this pitching moment to the nose of 
the decelerator (pin B) it is seen that 


S Z c l 

k, = * k 

Z S 7 s 
s s 


- 2 
it c 

S s “ IT - 


- 2 

q * C Z - 2 
s z = = * c s 


1^5 




where x is the transfer distance. In the nomenclature of the 
present investigation this is equivalent to a^. Thus, 



From equation (3 6 ) the spring constant is therefore 


k = 


- c C 


m a,B 


(38) 


It is obvious from this expression that for any body of revolution, 

C will be negative valued so that the spring constant will 

a,B 

always have a positive value. 

Natural angular frequency .- From equation (26) the natural 
angular frequency is expressed as 
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Assuming the decelerator to be a homogeneous solid cone, the terms 

m,-^ 2 + J 2 can s ^ own to be 

* 2*2 + J 2 = °2 ^ (if + k) (59) 

(see reference 9) 

Rewriting the expression for the natural angular frequency, 


CD 

n 


_ / a 2 N 

c ( — °N - C m 
\c a a> 


of -i £ 


(40) 


It can be seen from this expression that increasing the distance from 
the center of gravity to the apex of the cone (i.e. increasing a^) 
results in corresponding decreases in the natural angular frequency. 
Since a homogeneous distribution of mass for the cone has been 
assumed, a center of gravity shift would necessitate adding mass to 
the system. Equation (4-0) shows that increasing the mass twofold 
would cause a decrease in cd by a factor of i JW . 

It has been shown that if the diameter of the cone is doubled, 
the aerodynamic spring constant increased by a factor of 8. The effect 
of changing the cone diameter on the natural angular frequency will 
now be investigated. 

Let 



where 




xsarst 
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s s small cone 
Z = large cone 


Noting that & 2 / r wi H be the 8ame for any diameter of a cone having 
the same semiapex angle, the ii.ertia terns can be written as 

f 


T m „ 2 ( 16 *2 3 

‘s " Vs \15 + 20 y 


T _ — y. 2 ( 16 a 2 ,3 

X Z " Vl \15 2 + 20 

' r 


Assume the con° diameter to be doubled so that r 7 = 2 r j also 

4 S 


assume m, - 2 m . Thus, 
Z s 


+ ' 8I 


Therefore, 


CD = 

n z 


8 k 

£ 

WT. 


= CD 


s s 

It is seen that, if the mass of the cone is assumed to double along 
with the diameter of the cone, then there is no change in the natural 
angular frequency. On the other hand, if the mass did not double 
but was assumed to remain constant, then 


co = J2 CD 
n Z ' n s 

Another cone property, which might be considered along with the 


mass, is volume. Increasing cone diameter necessarily increases cone 
volume, so that it might be useful to consider two cones having the 
same density (p). 


40 


Density is expressed as 


m n, 

„ s l 

e \ v i 

1 2 

where cone volume, V, is given by j n r h. Writing 


.. 1 2 . 
V ! “ J” r i h I 


and considering the diameter and height of the cone to be doubled, it 
is seen that 


r = 2 r and h, = 2 h . 
Is Is 


Thus, 


8 ? 

V, = - it r h =8 V 
13 s s s 


Because of the condition of constant density, the mass of the larger 
cone is seen to be 


m, = 8 m . 

I s 

This leads to I 7 = I , so that the natural angular frequency of 

L S 

the larger cone is 





1 

2 


CD 

n 

s 


Stability boundary .- A final comment to be made about the motion 
of the decelerator about its nose concerns the stability. It has 
been stated that a criterion for stability is that the damping be 
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positive values. From the governing equation of motion, equation (24), 
it is seen that to insure stability it is necessary thet 

2 

a 2 3 l a 2 dM 3 m s 
V 3cc " V 3a ” 3q 

As mentioned previously for conical decelerator shapes, lift-curve 
slope is the aerodynamic parameter effected most by changing cone 
semiapex angle (see Figure 7) so that the boundary between stability 
and instability is given by 


3L 1 cM V °° 3M 
3cl = aT 5a + 5 3q 


(41) 


Motion of Rod and Decelerator About Pivot Point - 
Special Case Mo. 2 
For this particular case the A 

motion of the rod and decelerator about 
the pivot point (pin A) is considered, 
the rod and decelerator being rigidly attached at B. To obtain the 
equation of motion of 9 it is beneficial to refer back to the 
coordinate system shown in Figure 3 and to set up a free-body 
diagram as follows: 
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hi^h 


Thus, 


a = 


0 + 



(*3) 


Substituting a and a, into the expressions for the aerodynamic 
lift and pitching moment yields results similar to equations ( 9 ) and 
(10). 


= fl 4. i b A * a g) a ^ 5l ( b i + a s) 3 
5a 55 v_ H + rr 0 + r7 v 9 

oa oa o# 


m 


5m 


( f + » 2 1 

+ i - — 


m - ^ e + 1. - — ^ 0 + ^9 + ^0 + ^M§ 


b, + a, 


?) 5m a 


5a 


00 5a 5§ 


(*+5) 


where 


c3M _ 5 m 

as ’ *5 


Substituting these expressions into equation (42) and collecting terms 
yields the following equation of motion for the one degree of freedom 


system: 


+ 8 2)%L ( b l + 

> 2 ) 

5m 

^°° 5<x v » 


5a 

) & . + * 2 ) 

5m 

5m 

21 5a V . 

5a 

S& 




0 - m L ga 1 + mgg^ + a 2 ) (46) 


At- previously, terms involving ^ and — are neglected and this 

da dd 

equation reduces to 

- „ 2 A _ /v . . . , ] s , T( b l + a 2) dL ( b l + **2) dM SmI A 

Vl m 2l b l + a 2j + 0 + [_ iT 3£ r 3S " 5qJ 6 


(h + a 2 ) ( E + H) - 35 6 - "l 8 a l + “28 (\ + a 2 ) (47) 


The particular solution is found from this equation to be 


0 ^ 
P 


B 1 ga 1 + “gS (\ + a g ) 

l b l + a 2j ( D + 35 ) - 


This can be simplified somewhat by recalling that 


dL dN A . A o 

*r — — 'r— - A at a = 0 
da oa 


(48) 


Thus, 


”l ea l + ”2 8 (\ + a 2 ) 

Vir^rsTST 

( b l + a 2j 3S “ 3S 


The homogeneous equation is given by 
2 


2 r \2 

n^a^ + m 2 IK + a 2 J + J 


(*9) 


« f Ql * * 2 ) dL ( b l + a 2) dM dM 

Vco 3a V» 3a ~ 5q 


n 


( b i + * 2 ) ( D + H)- si 


0 = 0 


(50) 


which has the same form as equations (22) and (23), where 




2 

r 

Vi‘ 

’ + m s(\ + * 2 ] 


In 2 

[(5 + “ 2 ) 

1 dL ( b l + “ 2 ) dM 

Sa V Sa 

00 

dM 
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Spring constant .- From equation (50) the spring constant can be 
written as 


k = ( b i + *2) ( D + si) “ si 


Nondi mensional i zing this expression, 
k - <?»S c 




'L ~ da 
a 


Sc 


m 

m da 
a 


This can be written as 


k = 


9cc S c 


( b i ♦ *2) 


s - c 


m 

a a 


(52) 


The spring constant is seen to be positive valued at all times, the 
magnitude of which remains invariant with center of gravity location. 
Similar to the argument presented in describing the motion of the 
decelerator about its nose, equation ( 58 ), the spring constant in the 
above expression can be written as 


(53) 

a, A 

where is the aerodynamic pitching moment referenced to the 

a, A 

point of rotation, pin A. The magnitude of C (and thus the 

m_ . 



in equation 


spring constant) in this expression is larger than C 

m a,B 

(38) by the amount - — CL. . 

c a 

To evaluate the spring constant for two conical shaped decelera- 
tors having the same semiapex angle, the spring constant is written as 
follows: 
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where C, T and C are the same for the two cones. Assume 
N m 

a a 

b x +a 2 


is equivalent to 
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so that 
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which is the same expression obtained in equation (37) • Considering 
the larger cone to have twice the diameter of the small cone, the 
comments made concerning equation (37) are appropriate here. However, 
there is a condition placed on the length of the rod by the expression 
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Since 


then 


b l + a s X + a 2 

s s l l 
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Therefore,, when the diameter of the cone is doubled ^c ^ = 2 c 
the length of the rod for the larger cone must he twice that for the 
smaller cone ^b^ = 2 b^ j in order that equation ( 5 M be valid. 

Natural angular frequency . - From equation (50) the natural 
angular frequency is expressed as 
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(55) 


Assuming the rod and cone to be homogeneous solids, the mass moment 
of inertia (j) about their centers of gravity can be written as 
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where 
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and 
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Letting = m^, 
4 

h = ^ & 2 , and a 1 
expressed as 


which isn't unrealistic for a decelerator system, 
^1 

~ 2 ~ > the natural angular frequency can be 
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From this it is seen that increasing dynamic pressure causes an 
increase in o)^, while increasing mass and/or the distance of mass 
from the rotation point (pin A) causes corresponding decreases in 
o> n . Assuming the length of the rod to be much greater than the 
distance from the cone CG to the nose of the cone ^b^ > > a^^ 
will not seriously alter these trends so that equation (56) may be 
reduced to 



( 57 ) 


The effect of doubling the cone diameter can be seen by letting 



where 
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1 

It is noted that letting — be the same for the two cones leads to 
the same remarks made in the previous section (see page 47) concerning 
the effects of doubling the cone diameter on the natural angular 
frequency. With r ? = 2 r , it is seen that 




or b^ 



The reader will note that this is the same condition concerning the 
length of the rod which resulted during the previous discussion of 
the effects of doubling cone diameter on the spring constant. 

Stability boundary .- Since the spring constant is always positive, 
it is only necessary for the damping to be positive to insure stability 
for the dynamic system described by equation (47). Thus, the boundary 
between stability and instability is represented by 

( b l + a 2) dL ( b l + a 2) dM dM „ 

~^T 5S r — - 55 = 0 

0O 09 ^ 


( 58 ) 
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In terms of lift-curve slope this becomes. 


&L _ 1 dM V « bU 

33 = ( b l + a a) 35 (b x + a 2 ) 2 35 


(59) 


The similarity between this stability boundary and that obtained for 
the cone rotating about its nose, equation (4l), should be noted. 


VII. CONCLUDING REMARKS 


A preliminary investigation has been undertaken to theoretically 
determine the geometric and aerodynamic parameters which have the 
greatest Influence on the stability and performance of the nonporous 
towed decelerator. To aid in this endeavor a mathematical model was 
generated to describe the flexible tow line - decelerator dynamical 
system as a rigid two-body problem. The resulting second-order 
governing differential equations of motion were used to obtain the 
characteristic equation (quartic) describing the coupled motions of 
the two bodies. Evaluation of the coefficients of the quartic yielded 
expressions which illustrate when the dynamical system is unquestion- 
ably unstable. 

The two single-degree of freedom cases investigated were the 
motion of the decelerator about its nose and the motion of the rod 
and decelerator about the pivot point. For both of these cases 
expressions were derived for the natural angular frequency, damping 
factor, steady-state solution, and time to damp to one-half amplitude. 
Specific comments were made concerning a conical shaped decelerator 
to illustrate some of the geometric and aerodynamic parameters 
affecting the spring constant and the natural angular frequency. The 
boundary between stability and instability was obtained for these 
two single-degree of freedom cases. 
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